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We review the particle physics ingredients affecting the normalization, shape, 
^sO and flavor composition of astrophysical neutrinos fluxes, such as different pro- 

duction modes, magnetic field effects on the secondaries (muons, pions, kaons), 
__ and flavor mixing, where we focus on ^7 interactions. We also discuss the in- 

[T'j terplay with neutrino propagation and detection, including the possibility to 

detect flavor and its application in particle physics, and the use of the Glashow 
resonance to discriminate pry from pp interactions in the source. We illustrate 
the implications on fluxes and flavor composition with two different models: 
1) the target photon spectrum is dominated by synchrotron emission of co- 
accelerated electrons and 2) the target photon spectrum follows the observed 
+3 photon spectrum of gamma-ray bursts. In the latter case, the multi-messenger 

extrapolation from the gamma-ray fluence to the expected neutrino flux is 
i— i highlighted. 

> 

1 Introduction 

\o 

In addition to gamma-ray and cosmic ray instruments, neutrino telescopes, 
. such as IceCube [l] or ANTARES [2] , provide interesting data on the sources of 

the highest-energetic particles found in the universe, so-called "cosmic acceler- 
ators" ; see Refs. (3j|6] for reviews. In particular, neutrinos are a prominent way 
I to search for the origin of the cosmic rays, or to discriminate between leptonic 

• • and hadronic models describing the observed spectral energy distribution of 

. £^ photons. There are numerous possible sources, see Ref. [7] for an overview 

and Ref. [8] for the general theory. Very interesting extragalactic candidates 
for neutrino and cosmic ray production may be gamma-ray bursts (GRBs) |9 



and active galactic nuclei (AGNs) 10 12 . The most stringent bounds for 



these sources, which are expected to be roughly uniformly distributed over 
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the sky, so far come from IceCube, which has recently released data on time- 
point source searches, GRB neutrino 



integrated 13 and time-dependent 14 
searches [15], and diffuse flux searches 
have been detected, which has been for a long time consistent with generic 
Waxman-Bahcall 



16 . So far, no astrophysical neutrinos 



17 



18 bounds. How- 



and Mannheim-Protheroe-Rachen 
ever, data from IC40 and IC59, referring to the 40 and 59 string configura- 
tion of IceCube, respectively, start to significantly exceed these bounds, see 
Refs. 15 16 19 , which is in tension with the corresponding neutrino produc- 
tion models, such as Refs. [9 p0|2l] for GRBs. For example, neutrino data may 
soon challenge the paradigm that GRB fireballs are the sources of the ultra- 
high energy cosmic rays (UHECR) [22]. For constraints to AGN models, see, 
e.g., Ref. 23 . As a consequence, the age of truth has come for neutrino astro- 



physics, which is also the age of precision: Especially since data are available 
now, it is necessary to critically review the underlying assumptions from both 
the astrophysics and particle physics perspectives, and to develop the models 
from rough analytical estimates into more accurate numerical predictions. 

In this review, we focus on the minimal set of particle physics ingredients 
for the neutrino production, which must be present in virtually all sources, 
using several specific examples. We focus on photohadronic (jrf) interactions 
for the meson production, with the exception of Sec. |4.3| We do not only dis- 
cuss the predicted neutrino flux, but also the flavor and neutrino-antineutrino 
composition at source and detector. The discussed effects include: 

• Additional pion production modes, such as t-channel (direct) and multi- 
pion production; see, e.g., Refs. 24-30 . 



2!) 



36 



2 



Neutron and kaon production; see, e.g., Refs. 

The cooling and decay of secondaries (pions, muons, and kaons); see, e.g., 
Refs. " 



29 30 33 34 37 39 



Flavor mixing and possible new physics effects; see Ref. [40] for a review. 
The helicity-dependence of the muon decays; see, e.g., Refs. 28 38 ^] 
Spectral effects, such as the energy dependence of the mean free path of 
the protons, and their impact on the prediction; see, e.t 



Refs. 44 45 



The impact of the maximal proton energy on the neutrino spectrum; see, 
e.g., Ref. 36 . 

Deviations from the frequently used E~ 2 neutrino flux assumption; see, 
e.g., Ref. [46]. 



While many of these effects have been studied elsewhere in the literature, we 
mainly show examples generated with the NeuCosmA ( "Neutrinos from Cos- 
mic Accelerators" ) software in this review to present them in a self-consistent 
way. 

2 In general, there is also an additional contribution from charmed meson produc- 



tion, see Refs. [32133 



for a detailed comparison. 



3 This effect has been discussed earlier in the context of atmospheric neutrinos, see, 
for instance, Refs. 41-43 
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The structure of this review is as follows: In Sec. [2] we give a simplified pic- 
ture for the connection among neutrinos, cosmic rays, and gamma-rays. Then 
in Sec. |3j we review the minimal set of ingredients for neutrino production 
from the particle physics perspective. In Sec. |4j we discuss neutrino propa- 
gation and detection, including the possibility to detect flavor and the use of 
the Glashow resonance, where we illustrate how new physics can be tested 
in the neutrino propagation in Sec. |4.4| We furthermore present two specific 
applications: a generic AGN-like model in Sec. [5] and a model for GRBs in 
Sec. |6j where the main difference is the model for the target photons. Then 
we finally summarize in Sec. [7] 



2 Neutrinos and the multi-messenger connection 

Here we outline a simplified picture of the neutrino or cosmic ray source, 
as often used in the literature, whereas we add extra ingredients in the next 
section. In this approach, charged mesons originate from pp or pj interactions, 
where we focus on jry (photohadronic) interactions in this work; see, e.g., 



Refs. 47 48 for pp interactions, which may be dominant for particular source 
classes, such as supernova remnants. In the simplest possible picture, charged 
pions are produced by the Z\(1232)-resonance 

{n + 7T + | of all cases 
2 f u ■ W 

p + 7r u | of all cases 

While this process is not sufficient for state-of-the-art models for neutrino 
production, it is very useful to illustrate a few qualitative points common to 
many cosmic ray and neutrino production models. The protons on the l.h.s. of 
Eq. ([I]) are typically assumed to be injected into the interaction volume with 
an {E' p )~ a spectrurrpl coming from Fermi shock acceleration, where a ~ 2. 
They interact with the photons on the l.h.s. of Eq. ([I]) with energy e' ~ 
(0.2 — 0.3) GeV/E' p . While the assumptions for the injected protons are similar 
for most models (except from the minimal and maximal energies), the target 
photons are typically described in a model- and source-dependent way, such 
as by: 

1. synchrotron emission from co-accelerated electrons or positrons, 

2. thermal emission, such as from an accretion disk, 

3. a more complicated combination of radiation processes, 

4. an estimate inferred from the gamma-ray observation, 

just to name a few examples. While any realistic simulation of a particular 
source will imply option 3) , the other options typically rely on fewer parame- 
ters and may be good approximations in many cases. In particular, a reliable 



4 Here primed parameters refer to the shock rest frame (SRF), where as unprimed 
parameters to the observer's frame. 
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prediction for the photon density in the source may be obtained from the 
gamma-ray observation, option 4), if the photons can escape. In fact, we will 
use option 1) in Sec. [5] and option 4) in Sec. [6] 

After an interaction between proton and photon, the particles on the r.h.s. 
of Eq. are produced with the given branching ratios. The neutrinos then 
originate from ir + decays via the decay chain 

7T+ -> /!+ + Ufj, , 

/i+ -> e + + v e + , (2) 

where in this standard picture v e : 17, : v T are produced in the ratio 1:2:0 
if the polarities (neutrinos and antineutrinos) are added. In addition, high- 
energy gamma-rays are produced by 

98.8% 



7T * 



7 + 7. (3) 



These are typically emitted from the source at lower energies due to electro- 
magnetic cascades, in addition to gamma-rays escaping from the interaction 
volume (the ones contributing on the l.h.s. of Eq. Q). 

From Eq. 0, we can also illustrate the production of cosmic ray protons, 
ignoring for the moment that the composition of cosmic rays may be heavier at 
high energies [49] . First of all, some of the protons injected into the interaction 
volume on the l.h.s. of Eq. may escape, leading to cosmic ray production. 
However, even if the protons are magnetically confined, the neutrons on the 
r.h.s. of Eq. 0, which are electrically neutral, can easily escape if the source 
is optically thin to neutron escape. After decay (typically outside the source) 

n — > p + eT + D e , (4) 

they lead to cosmic ray flux and an additional D e neutrino flux which is an un- 
avoidable consequence of the interactions in Eq. 0. The cosmic ray protons 
with energies above 6 • 10 19 eV interact with the cosmic microwave background 
(CMB) photons by Eq. 0, leading to the so-called Greisen-Zatsepin-Kuzmin 
(GZK) cutoff [50,51 . However, according to Eq. 0, charged pions are pro- 



duced in these interactions as well, which means that an additional neutrino 
flux should come with that, which is often called "cosmogenic neutrino flux". 

In summary, the photohadronic interaction in Eq. offers a self-consistent 
picture for a cosmic ray source, with a possible connection among cosmic ray, 
neutrino, and gamma-ray escape. In specific models, however, that does not 
mean that a large neutrino flux is guaranteed for every cosmic ray source. For 
instance, the interaction rate for the process in Eq. 0, which depends on the 
photon density, may be low. 



3 Simulation of neutrino sources 



Here we give a more detailed generic picture of the simulation of neutrino 
sources from the particle physics perspective with the minimal set of ingredi- 
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ents, without using a specific model. A flowchart summarizing the contents of 
this section, which can be followed during the reading, is given in Fig. [6] 

3.1 Photohadronic interactions 

In order to describe the processes within an interaction volume (one zone 
in the simplest case), two kinds of spectra are needed: Q'(E') (in units of 
GeV -1 cm~ 3 s _1 ) describes the number of particles injected or ejected per 
volume element and energy interval, and N'(E') (in units of GeV _1 cm -3 ) 
describes the particle density per energy interval. The secondary meson in- 
jection rate Qb(-Eft) for a pion or kaon species b produced in photohadronic 
interactions is given by (following Ref. [28] ) 

oo oo 

/dE' r 
-^N' p (E' p ) J cds'%(e')R b (x,y). (5) 

E ' b p o 

Here x — E' b /E' p is the fraction of energy going into the secondary, y = 
(E' p e') /mpQand Rb(x, y) is the "response function". If many interaction types 
are considered, the response function can be quite complicated. However, if it 
is known from particle physics, Eq. ^ can be used to compute the secondary 
injection for arbitrary proton and photon spectra. The important point here 
is that the secondary production depends on the product normalization of 
the proton density N' p (E' p ) and the target photon density A^(e') within the 
interaction volume. Thus a higher proton density can be compensated by a 
lower photon density, and vice versaj^] Another implication of Eq. ^ is that 
the secondary production depends on the densities within the source N' , not 
the injection rates Q' . Of course, one cannot look into the source, but can only 
observe cosmic messengers escaping from the source. As we will demonstrate 
later, the observed/ejected photon or cosmic ray spectrum Q' is only directly 
representative for the corresponding density spectrum within the source N' 
if "trivial" escape is the leading process, i.e., Q' = N' /t' csc with t' csc ~ R jc 
and R! the size of the interaction region. For this section, Eq. ([5| is used 
as a starting point for the computation of the neutrino fluxes, where we do 
not discuss the origin of spectral shape and normalization of N' and N'. In 
practice, typically an (E' p )~ 2 injection spectrum is assumed for the protons, 
as mentioned above, where the maximal energy is limited by synchrotron and 
adiabatic losses. The photon density may be a consequence of a complicated 

5 Here y can be related to the center-of-mass energy by s = m,p + 2m p (l — cos # P7 ) y, 
where 6 P1 is the angle between proton and photon momentum; 6 V1 — tt corre- 
sponds to heads-on collisions. 

6 Strictly speaking, this degeneracy only holds in the absence of any other radiation 
process. For example, inverse Compton scattering depends on the photon density 
individually, which means that the input spectral shapes (especially Ni, in Eq. |5jl) 
will be modified if that process contributes significantly. 
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Fig. 1. Contributions of different n production modes to the total spectrum for a 
typical GRB, an AGN, and a (10 eV) black body (BB) target photon density. Figure 
taken from Ref. [28]. 



interplay of radiation processes. In either case, the derivation of these densities 
depends on the model, and we will show several examples in Sees. [5] and [6] 

Compared to a numerical approach, Eq. ([I]) has limitations to describe 
the meson production. First of all, additional pion production modes con- 
tribute, such as higher resonances, direct (t-channel), and multi-pion produc- 
tion, which will also lead to 7r~ production (c/., Fig. [6]). These are not as easy 
to describe as the Z\-resonance because of different shapes of the cross sections 
and more complicated kinematics. The Monte Carlo package SOPHIA 25 can 



deal with these interactions. In order to increase the efficiency, often parame- 



terizations of SOPHIA are used, such as Refs. 28 52 . In the following, we use 
Ref. [28] (model Sim-B), because the secondary muons and pions are needed 
explicitly. We show the impact of the resonant production (including higher 
resonances) on ir + production in Fig. [I] for a typical GRB, an AGN, and a 
(10 eV) black body (BB) target photon field. As one can read off from this fig- 
ure, the resonances always give a reasonable first estimate for the actual pion 
production, but quantitatively they only dominate at the breaks. In addition, 
multi-pion processes can change the spectral shape significantly, such as for 
the GRB example, which is a consequence of the cross section dependence 
on the center-of-mass energy. As a further limitation, note that Eq. ([I]) does 
not describe kaon production and subsequent decay into neutrinos, where the 
leading modes are given by 

P + -f^K + +A/S, (6) 
K+ -> /!+ + i/„ . (7) 

The branching ratio for the leading channel in Eq. Q is about 64%. The 
second-most-important decay mode is if ± — > ^ + n" (20.7%). The other 
decay modes account for 16%, no more than about 5% each. Because interest- 
ing effects can only be expected in the energy range with the most energetic 
neutrinos, we only use the direct decays from the leading mode. 
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In the literature, the Zl-resonance approximation in Eq. ([I]) is, even in 
analytical approaches, typically not taken literally. For example, a simple case 
is the approximation by Waxman and Bahcall (9| ("WB Z\-approx."), for 
which one can write the response function as 



R„±(x,y) — 0.5x<5(x— 0.2)x500/ibarn 



x < 



2y< 0.2 GeV 

l-fo^v) 2 0.2GeV<2, 

< 0.4 GeV 



0.12 GeV 
(2y) 2 



2y > 0.4 GeV 



which implies that charged pions are produced in 50% of all cases, and these 
take 20% of the proton energy. In addition, the width of the Z\-resonance 
is taken into accountQln fact, this function peaks at 2y ~ 0.4 GeV, which 
is higher than the threshold for photohadronic interactions - and it is even 
a little bit higher in the numerical calculation. The reason is that for the 
threshold often head-on collisions are assumed {6 pl — tt), whereas these only 
contribute a small part to the total number of interactions. Using Eq. ^ in 
Eq. ([5]) and re-writing the integral over e' in one over y, it is easy to show 
that for power law spectra N' p oc (E' p )~ a and N' oc (e') ^ 



,-at+P-l 



(9) 



This means that the pion spectral index depends on both the proton and 
photon spectra, where it is inversely proportional to the photon spectral index. 
As we will see below, the neutrino spectrum follows the pion spectrum, which 
means that the assumption of an E~ 2 spectrum for the neutrinos, as it is 
often used in data analysis and many models in the literature, is only a valid 
assumption for j3 ~ 1 - which is roughly observed for GRBs below the break. 
On the other hand, if the target photons come from synchrotron emission, such 
a hard photon spectrum is not possible, and (3 > 3/2 may be more plausible 
if the electrons are injected with a spectral index similar to the protons. As a 
consequence, the neutrino spectrum becomes harder. In addition, multi-pion 
processes in the photohadronic interactions will act in the same direction and 
make the neutrino spectrum even harder, c/., Fig. [I] see also Ref. |29 for a 
detailed comparison between the approximation in Eq. pi and the numerics. 
Note that for pp interactions with "cold" (non-relativistic) protons, the E~ 2 
assumption may be plausible 47 . We discuss the implications of the E~ 2 



assumption for the detector response in Sec. |4j 



3.2 Decays of secondaries 

The weak decays of pions and muons are described in detail in Ref. [38] . In 
general, in case of ultra-relativistic parents of type a, the distribution of the 

7 Note that this description is slightly more accurate than Ref. [9], which uses an 
additional integral approximation. 
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Fig. 2. Flux ratio at the source between electron and muon neutrinos as a function 
of the neutrino energy, for the same objects as in Fig. [I] No secondary cooling is 
included in this figure. Figure taken from Ref. |28|. 



daughter particle of type b takes a scaling form in order to obtain for the 
energy spectra 

Q'b( E 'b) = E ^ dE 'a KiK) (Cc)" 1 ^7 ( j|) (10) 

summed over all parent species. The functions F a ^ for pion, kaon and hclicity 
dependent muon decays can be read off from Ref. [38] (Sec. IV). Consider the 
simplified case of a 8- function for F a ^b. For instance, for decays of neutrons 



in Eq. Q, one may approximate 
with Xn — 5.1 x 10 -4 . If the neutrons do not interact, Q' n — ) _1 ( see 



Xn) (ID 



below), and we find for the neutrino injection Eq. (10) 



Qk(K) = -Q'n(—) • (12) 

Xn \XnJ 

In this case, the neutrino spectrum follows the neutron spectrum. If the neutri- 
nos originate from pion decays, the neutrino injection follows the pion injection 
spectrum by similar arguments. 

Since the decays of muons are helicity dependent, i.e., F a ^b is different 
for left- and right-handed muons, it is necessary to keep track of these two 
species separately (c/., Fig. [6]). Although the effects of the helicity-dependent 
muon decays on the fluxes are probably small, the flavor composition is slightly 
affected, depending on the parameters of the source. We illustrate this effect in 
Fig. [2] for the GRB, AGN, and black-body examples in Fig. [T] In this figure, 
the horizontal line corresponds to the standard assumption, i.e., neutrinos 
being produced in the flavor composition v e : : v r of 1 : 2 : 0. Including 
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the scaling of the secondary decays, the dotted (blue) curves are obtained 
if the helicity of the muons is averaged over. From the comparison with the 
light (red) solid curves, it is clear that this assumption is implemented in 



SOPHIA. On the other hand, as pointed out in Ref. 38 , keeping track of 
the muon helicity slightly changes the flavor composition, see dark (black) 
solid curves, which are significantly different from the standard assumption 
and the helicity- averaged version. However, it is also clear from Fig. [2] that 
the deviation from the standard prediction depends on the input spectra, as 
it is smaller for the AGN than the GRB example. The dashed (green) curves 
show the contributions of the neutron and kaon decays, which affect the flavor 
composition at very low and high energies, respectively. 

3.3 Cooling of secondaries 

In order to describe the cooling of the secondary pions, muons, and kaons, 
we use the steady state approach, i.e., we do not allow for an explicit time- 
dependence since the statistics of neutrino observations is typically expected 
to be low. The steady state equation for the particle spectrum, assuming 
continuous energy losses, is given by 

d N'(E') 

Q ' {E ' ] = Je> {b ' {El) N {E )) + ' (13) 

with t' esc (E') the characteristic escape time, b'(E') = —E' t[~^ s with t[~^ s (E') = 
— 1/E' dE' jdt' the rate characterizing energy losses. This differential equation 
balances the particle injection on the l.h.s. with energy losses and escape on the 
r.h.s. of the equation. Note again that the steady density N'(E') is needed for 
the photohadronic interactions in Eq. (J5|, not the actual injection spectrum. In 
addition, note that if there are no energy losses (6' = 0), one has immediately 
Q'(E') = N'{E')/t' csc from Eq. ([13]), which we have already used above. If 
decay is the dominant escape mechanism, one finds Q'(E') cx N' (E')/E' . 

While the primary proton and photon spectra in Eq. |5]) could be affected 
by a number of radiation processes, for the neutrino fluxes and flavor com- 
positions, at least the processes of the secondaries (pions, muons, kaons) are 
important. We illustrate the synchrotron cooling and decay rates for pions, 
muons, and kaons as a function of energy in Fig. [3] As one can easily see in 
the figure, for any species, decay dominates at low energies, while synchrotron 
cooling dominates at high energies. Other cooling or escape processes are of- 
ten sub-dominant, as illustrated by the gray-shaded region for an adiabatic 
cooling components in GRBs. The two curves meet at a critical energy E' c 
for each species, which is different depending on the particle physics param- 



eters. As a consequence of Eq. (13), the corresponding steady spectra N' are 



loss-steepened by two powers above 



_ / 97re m 5 c 7 
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Fig. 3. Synchrotron cooling and decay rates for pions, muons, and kaons as a 
function of energy for a magnetic field of B' ~ 300 kG. The gray-shaded region 
shows the estimated range for the escape (or adiabatic cooling) rate for GRBs for a 
variability timescale ranging from 0.01 s to 1 s. Figure taken from Ref. 30 



where synchrotron and decay rates are equal. These critical energies depend 
on the particle physics properties of the parent, i.e., the mass m and the 
rest frame lifetime to, and the magnetic field B' as the only astrophysical 
parameter. It is therefore a very robust prediction, and might allow for the 
only direct measurement of B' . Re-scaling the magnetic field shifts the critical 
energies by a constant amount on the horizontal axis, but does not change 
the spacing in the logarithmic picture. In Fig. [3j we also show the estimated 
range for adiabatic cooling as additional cooling component (shaded region), 
which may have some impact, especially on the muons, in extreme cases. In 
these cases, the height of the spectral peaks will be somewhat reduced, but 
the qualitative picture does not change. 

We show in Fig. [4] the consequences for the spectral shape at one GRB 
example. Here the injection spectrum of electron neutrinos (left panel) and 
muon neutrinos (right panel) is shown, including the individual contributions 
from the parents. First of all, one can read off from this figure that in the 
case of no losses (dashed curves) the spectral shapes of all contributions are 
very similar, and the neutrino fluxes add in a trivial manner. A change of the 
primary spectra N' p and N' y in Eq. ([5| may change the shape of the dashed 
curves, but almost in the same way for all curves. If the synchrotron losses 
are switched on (solid curves), the spectral split predicted by Eq. (14) (see 
also Fig. [3]) among the neutrino spectra coming from different parent species 
can be clearly seen in the right panel. One can also see a small pile-up effect 
coming from the muon decays, i.e., a small region where the cooled muons 
coming from higher energies pile up and dominate, and lead to a higher flux 
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Fig. 4. Injection spectra for electron neutrinos (left panel) and muon neutrinos 
(right panel) at the source for a GRB example, where the individual contributions 
from the parents are shown (neutrinos and antineutrinos added). Figure taken from 
Ref. [30] . 



than in the "no losses" case. In the left panel, only two spectra are shown, since 
only muon or neutron decays may produce the electron flavor. Because of the 



very small Xn in Eq. (12 1, the neutron decays only show up at low energies. 
Since the decay of pions, muons, or kaons, which dominate in different energy 
ranges, lead to different neutrino flavor compositions, the flavor composition 
at the source will be changed as a function of energy. In the literature, often 
the following source classes in terms of u e : : v T are distinguished^] 

Pion beam source Pion decays and muon decays equally contribute, flavor 

composition 1:2:0. 
Muon damped source Strong muon cooling, which means that pion decays 

dominate; flavor composition 0:1:0. 
Muon beam source Pilc-up muons dominate; flavor composition 1:1:0. 
Neutron beam source Neutron decays dominate; flavor comp. 1:0:0. 
Undefined source Several of these processes compete; flavor composition 

X:1-X:0 (X£{0,i,i,l}). 

Note that in none of these cases v r are produced in significant quantities, 
which means that it is sufficient to use the ratio between the v e and pro- 
duction at the source, sometimes called "flavor ratio". From the preceding 
discussion, it is clear that the above classifications can only hold in specific 
energy ranges. However, all of these sources can be recovered in a numerical 
simulation. In Fig. [5j four examples for energy dependent flavor ratios at the 



See Ref. 53 for a summary, including other options leading to a neutron beam 



or muon beam-like source. 
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Fig. 5. Energy-dependent flavor ratio at the source for several selected parameter 
sets (test points "TP" ; c/., Fig. 10 1 for a model where is generated by synchrotron 
losses of co-accelerated electrons. The (dark) gray-shaded areas mark the regions 
where the fluxes peak. Figure taken from Ref. [34|. 



source are shown for different parameter sets, where in this model is gen- 
erated by synchrotron losses of co-accelerated electrons. In this figure, also 
the relevant flavor ratio ranges for the different sources introduced above are 
shown. The upper right panel shows the classical pion beam source, which is 
typically found for low magnetic fields. Nevertheless, the contribution of neu- 
tron decays at low energies can be clearly seen. The lower right panel shows 
a pion beam evolving in a muon damped source at high energies, as, e.g., in 



Ref. 37 . The upper left panel depicts a muon beam to muon damped source. 
In this case, the cooled muons pile up at lower energies, where the muon de- 
cays dominate. And the lower left panel shows an undefined source, where 
several processes compete. 

Of course, not only the secondaries are affected by synchrotron (or adi- 
abatic) losses, but also the primary protons. Depending on the model, one 
can use these losses to determine the maximal proton energy, or one can 
put in the maximal proton energy by hand. From Fig. [4j it is interesting to 
discuss the impact of the maximal proton energy on the neutrino fluxes. In 
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the "no losses" case (dashed curves), the maximal neutrino energy is directly 
determined by the maximal proton energy E' l/max ~ 0.05 E' max . There is, 
however, one exception: the neutrinos from neutron decays are limited by 
E' v ax ~ 10 _4 £'p max , c/., Eq. (12 1. If the synchrotron losses are switched 
on (solid curves), the neutrino spectrum from neutron decays still follows 
the proton spectrum, since the neutrons are electrically neutral, whereas the 
maximal neutrino energies for the other production modes are determined by 
Eq. (14). As a consequence, the neutron decay spectrum strongly depends on 
the assumptions for E' p mltx , whereas the other spectral shapes are entirely 
unaffected by E' max as long as E' max >6E' C for kaons (see Eq. (14 1, factor 
six from kaon production and decay kinematics). It is therefore not surprising 
that for strong enough magnetic fields one can find parameter sets for which 
the neutron decays dominate (c/., Refs. [35|[36] ), but one should keep in mind 
that this depends on the assumptions for the maximal proton energy (and the 
inclusion of multi-pion processes etc., which may mask this effect). 



3.4 Transformation into observer's frame 



The transformation of the injection spectrum of the neutrinos Q' Va from the 
source to the observable flux <j)p of vp (in units of GeV -1 cm -2 s _1 ) at the 
Earth is given by 

h= £ NP aP { ±±*lQ' Va , E v = ^-E' v , (15) 

a—e,fi,r L 

where a simple Lorentz boost r is used (instead of a viewing angle-dependent 
Doppler factor). Here P a p is the transition probability v a — > i>p, discussed in 
Sec.^ and TV is a (model-dependent) normalization factor. For example, if an 
isotropically emitting spherical zone is boosted with r towards the observer, 
then N = (4/3) R' 3 irr 2 since the emission is boosted into a cone with opening 
angle l/T. For a relativistically expanding fireball, it is simpler to perform 
the transformation in a different way, see Sec. [6j Furthermore, cIl(z) = (1 + 
z) d com (z) is the luminosity distance, and d com (z) is the comoving distance. 



From Eq. (15), one can read off that the rcdshift dependence of the neutrino 
luminosity scales as E^cf) cxl/d\ independent of the model, which is expected. 



Note that in Eq. (15) the neutrino and antineutrino fluxes are often added if 



the detector cannot distinguish these. 



3.5 Summary of ingredients, and limitations of the approach 

We summarize the generic neutrino production in Fig. [6] As one can see in 
this flowchart, the starting point are the proton and photon densities within 
the source. Once these (and B') are fixed, the rest is just a particle physics 
consequence. Therefore, for the computation of specific neutrino fluxes, the 
main effort is actually to determine N'(E' ), N'(e'), and B' . 
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Fig. 6. Flowchart describing the model (in SRF). The functions Q'(E) denote 
(injection) spectra per time frame [(GeVcm 3 s) 1 ] and N'(E) steady spectra 
[(GeVcm 3 ) 1 ] derived from the balance between injection and losses or escape. 
Dashed arrows stand for solving the steady state differential equation Eq. (13 1. 



Figure taken from Ref. 30 



Of course, there are some processes not taken into account in this picture, 
which may add to the ingredients discussed above for specific source classes. 
For instance, secondary neutrons, produced in the photohadronic interactions, 
may interact again if the source is optically thick to neutrons, the secondary 
pions, muons, and kaons may be re-accelerated [54], synchrotron photons of 
the secondaries may add to N', etc. In addition, the neutrino spectrum may 
be more complicated in multi-zone models, since these naturally allow for 
more freedom. 

From a particle physics perspective, additional kaon and charmed meson 
production modes may be added, and also the secondaries may interact again 



(see Refs. 32 33 ). However, compared to analytical computations, the nu- 
merical approach described in this section already takes into account the sec- 
ondary cooling in a self-consistent way, additional neutrino production modes 
can be easily included, and the full energy dependencies can be accounted 
for. For example, it has been demonstrated in Ref. [44], that all necessary 
ingredients to reproduce the analytical GRB fireball (neutrino) calculations 
m Refs. [9j[l5j[20j[2l] are contained. It should represent the minimal set of 
of ingredients for neutrino production which are present in every source in 
the spirit of constructing the simplest possible model first. Of course, if B' is 
small, the secondary cooling effects will be small as well, which is automati- 
cally included. 
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4 Neutrino propagation and detection 

In this section, we discuss several aspects of neutrino propagation and detec- 
tion from the theoretical perspective. 



4.1 Neutrino propagation and observables 

It is well known that neutrinos may change flavor from the production to the 
detection point. While this phenomenon is in general described by neutrino 
oscillations, astrophysical neutrinos are typically assumed to suffer from deco- 
herence over very long distances. This means that effectively (in most practical 
cases) only flavor mixing enters the astrophysical neutrino propagation (see 
Ref. 55 for a more detailed discussion). In that case, P a p in Eq. ( 15 ) becomes 



\u a .. 



(16) 



for three active neutrinos, where U a % are the usual PMNS mixing matrix ele- 
ments in the standard parameterization; see e.g. 



Ref. 56 



for recent values of 

the mixing angles. This implies that neutrino oscillations, i.e., the Am 2 L/E- 
dependence, are averaged out. An initial flavor composition v e : : v T of 
1:2:0 will therefore evolve (approximately) into 1 : 1 : 1 at the detector, 
Refs. 



e.g., 



57 58 . In Sec. |4.4| we will see that Eq. ( 16 ) can significantly 



change in the presence of new physics effects, which opens new possibilities 



to test such effects. However, Eq. ( 16 ) also implies that there could be some 



sensitivity to standard flavor mixing, which may be complementary to Earth- 
based experiments, see discussions in Refs. 59-71 . In the light of the current 



bounds for astrophysical neutrino fluxes from IceCube, however, such appli- 
cations might be unlikely. 

The main observable in neutrino telescopes are muon tracks from charged 
current interactions of muon neutrinos, producing Cherenkov light, which can 
be detected in so-called digital optical modules (DOMs). Because of the long 
muon range which is increasing with energy, the muon track does not have 
to be fully contained in the detector volume, which leads to excellent statis- 
tics increasing with energy. In addition, muon tracks have a very good di- 
rectional resolution (order one degree). Additional event topologies include 
electromagnetic (mostly from electron neutrinos) and hadronic (from tau neu- 
trinos) cascades, as well as neutral current cascades for all flavors. For even 
higher energies, the tau track may be separated, leading to so-called double 
bang or lollipop events; see Ref. 



72 



for an overview. In practice, the main 
"flavor" analysis so far performed by the IceCube collaboration has been a 
cascade analysis [73]. To see that, consider that electromagnetic (from v e ) and 
hadronic (from v T ) cascades cannot be distinguished. A useful observable is 
therefore the ratio of muon tracks to cascades 59 
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Fig. 7. Flavor ratio R as a function of E v at the detector for a GRB neutrino flux 
with the Waxman-Bahcall (WB) "plateau" (in E 2 ^) between 10 5 and 10 7 GeV (c/., 
Ref. 17 ). The shaded regions show the impact of the 3cr mixing angle uncertainties 
, in about 2015 (next generation, dominated by Daya Bay and 
81 ), and in about 2025 (for a neutrino factory [82||83| ). Here #13 — is 
assumed for the sake of simplicity. Figure taken from Ref. 1 30 . 



now ("current 
T2K " 



R 



(17) 



Note that neutral current events will also produce cascades, and v T will also 
produce muon tracks in 17% of all cases, which, in practice, have to be included 
as backgrounds. In Ref. 



73 



cascade rate for a E„ 



the contribution of the different flavors to the 
extragalactic test flux with equal contributions of 
all flavors at the Earth was given as: electron neutrinos 40%, tau neutrinos 
45%, and muon neutrinos 15% (after all cuts). This implies that charged 
current showers dominate and that electron and tau neutrinos are detected 
with comparable efficiencies, i.e., that Eq. (17) is a good first approximation 



to discuss flavor at a neutrino telescope. The benefit of this flavor ratio is 
that the normalization of the source drops out. In addition, it represents the 
experimental flavor measurement with the simplest possible assumptions. For 
a pion beam source, one finds R ~ 0.5 at the detector, for a muon damped 
source R ~ 0.6, and for a neutron beam source R ~ 0.3, with some dependence 
on the mixing angles. In principle, this and other observables, such as the 
ratio between electromagnetic and hadronic showers, can therefore be used to 



determine the flavor composition of the sources, see, e.g., Refs. 74 -78] (unless 
there are new physics effects present, see e.g. Ref. 79 , which may yield similar 



ratios at the detector for different flavor compositions at the source). 

We show the ratio R as a function of energy for a GRB neutrino flux 
111 Fig. [7) At the "plat eau" of the flux, where the largest number of events 
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is expected from, one can clearly see the flavor transition between a pion 
beam and a muon damped source at higher energies. In addition, the neutron 
decays have an impact at very low energies. In this figure, the uncertainty 
on R (3cr) coming from the current mixing angle uncertainties is shown, as 
well as the expectation for the next generation of reactor and long-baseline 
experiments ("2015"), and for a high-precision neutrino oscillation facility 
("2025"). Obviously, the current uncertainties on the mixing angles are still 
too large to allow for a clear identification of the flavor composition at the 
source, whereas already the knowledge from the next generation will allow for 
a flavor ratio discrimination - at least in principle. 



4.2 Detector response and impact of spectral shape 



For time-integrated point source searches in IceCube 13 , the simplest possible 



approach to describe the event rate of muon tracks is to use the exposure 
Exp(E v ,6) = A° s (E„,5)t e x P , where A° ff is the neutrino effective area and 
t cxp is the observation (exposure) time. Here A^(E V ,5) is a function of the 
flavor or interaction type (which we do not show explicitely) , the incident 
neutrino energy E, and the declination of the source 6. The neutrino effective 
area already includes Earth attenuation effects (above PeV energies) and event 
selection cuts to reduce the backgrounds, which depend on the type of source 
considered, the declination, and the assumptions for the input neutrino flux, 
such as the spectral shape. Normally, the cuts are optimized for an E~ 2 flux. 
The total event rate of a neutrino telescope can be obtained by folding the 
input neutrino flux with the exposure as 



N = 



J dE l/ Exp(E 1/ ,6)<t> fi (E u ) = J dE„Af(E Ul S)t cxp ^{E u ). (18) 



Here (fi^E^) is, for point sources, given in units of GeV -1 cm~ 2 s _1 for neutri- 
nos and antineutrinos added. If backgrounds are negligible, the 90% (Feldman- 
Cousins) sensitivity limit Kqq for an arbitrarily normalized input flux used in 
Eq. (18) can be estimated as K go ~ 2.44/iV [84]. This imples that a predicted 



flux at the level of the sensitivity limit, irrespective of the spectral shape, 
would lead to the same number (2.44) of events. The 90% confidence level 
differential limit in terms of Efy^ can be defined as 2.3 i? l/ /Exp(£' ;y , S), see, 



e.g., Ref. 85 



The interplay between spectral shape and detector response has, for in- 
stance, been studied in Ref. [46) . In order to illustrate that, we show in Fig. [8] 
the limits for selected muon neutrino spectra and two different source declina- 
tions, corresponding to downgoing and upgoing events, for the time-integrated 
point source search in IC40. In this figure, the thick horizontal lines show the 
limits for an E~ 2 flux (in the dominant energy range), and the thick curves 
the differential limits. One can easily see that the differential limits are use- 
ful, since any "well behaved" neutrino flux will stay below these limits. In 
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8. Limits for selected muon neutrino spectra (including flavor mixing) from 



34 for differen t de clination bands for the time-integrated point source search 

2 flux (in the 



in IC40 (90% CL) [13]. The thick lines show the limits for an E, 
dominant energy range), and the thick curves the differential limits. Figure taken 
from Ref. [46]. 



addition, the differential limit shows the energy range where the instrument 
is most sensitive to a flux, whereas for the horizontal lines, representing an 
E~ 2 flux limit, only the contribution close to the differential limit minimum 
contributes to Eq. (18 1. In Fig. [8] the minimal and maximal energies for the 



horizontal lines are indeed arbitrarily chosen, since the limit hardly depends on 
thesej^From the different fluxes in Fig. [8] it is clear that the interplay between 
spectral shape and detector response is important. For example, spectrum #4 
will be better constrained in the left panel (downgoing events) than in the 
right panel (upgoing events) because the differential limit peaks at higher en- 
ergies - in spite of the lower absolute performance at the differential limit 
minimum. The reason is the coincidence between spectral peak and differen- 
tial limit minimum. This picture changes completely in the presence of strong 
magnetic field effects on the secondaries, see spectrum #2, where even the 
imprint of these effects becomes important. The comparison to the E~ 2 test 
flux (horizontal line) clearly demonstrates that the E~ 2 assumption is insuffi- 
cient for all practical purposes, since it does not take into account the energy 
range of the flux. In addition, as we discussed around Eq. |9| earlier, the E~ 2 
assumption only holds for a very special case for the target photon density. 
Therefore, detector response and source model are intimately connected, and 
it is not quite clear if some sources may be missed just because the detector 
response does not match the source prediction. We will come back to this 
issue in Sec. [5] 

9 For the chosen example, this holds as long as the energy range around the differ- 
ential limit minimum within about 1.5 orders of magnitude in energy is included. 
Sometimes the range where 90% of the events come from is shown. 
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Fig.[S]is also useful to illustrate the impact of the source declination, which 
is, for IceCube, a measure for the direction of the muon track in terms of the 
nadir angle because of the location at the South Pole. Note, however, that 
for other neutrino telescopes, such as ANTARES, this relationship is non- 
trivial. Obviously, the differential limit is different in the left and right panels 
of this figure, corresponding to downgoing and upgoing events, respectively. 
First of all, note that the downgoing events have to fight the atmospheric 
muon background, which leads to a worse performance at low energies be- 
cause of appropriate cuts. However, the upgoing events suffer from the Earth 
attenuation for energies PeV, which leads to the better sensitivity of the 
downgoing events for high energies, and the shift of the differential limit min- 
imum to higher energies. From the discussion above it is clear that both event 
types (upgoing and downgoing) are complementary not only because they test 
a different part of the sky, but also because they test different energy ranges. 

Finally, let us briefly comment on the statistics expected for neutrinos. 
First of all, it is clear from the discussion above that per definition any cur- 
rent IC40 limit, obtained over about one year of data taking, is compatible 
with about 2.4 events at the current 90% confidence limit. Assuming that 
the effective area increases by about a factor of three to four from IC40 to 



IC86 86 , one can extrapolate that the current limits are compatible with 
2.4 -4-10 ~ 100 events over ten years of full IceCube operation if the cur- 
rent bounds are saturated. Any further non-observation of events will reduce 
this maximal expectation. Therefore, one can already say right now that any 
conclusion about the astrophysical neutrino sources, the sources of the cosmic 
rays, or leptonic versus hadronic models for 7-ray observations will most likely 
be based on the information from many sources of one class. A typical exam- 
ple is the stacking of GRBs using their gamma-ray counterparts, such as in 



Ref. 15 . The aggregation of fluxes, no matter if diffuse or stacked, will how- 



ever imply new systematics and model-dependent ingredients, see discussion 



in Ref. 30 



4.3 Glashow resonance to discriminate pp from P7? 

A useful observable may be the Glashow resonance & e +e~ — > W~ — > anything 
at around 6.3 PeV 34 57 ,87 91 to distinguish between neutrinos and antineu- 



trinos in the detector, since this process is only sensitive to v e . For photo- 
hadronic (257) interactions, however, mostly 7r + and therefore v e are produced 
at the source, see Eq. ([I]) and Eq. pi), which means that no excess of events 
should be seen at this resonance energy - at least in the absence of flavor 
mixing. On the other hand, for pp interactions in the source, 7r + and tt~ 
are produced in about equal ratios, which increases the u e production rate. 
Therefore, the Glashow resonance is frequently proposed as a discriminator 
between pp and pry sources. Note that this argument is especially interesting 
for E~ 2 neutrino fluxes, whereas any other, (significantly different) spectral 
index may be a clear indicator for pr/ interactions. 
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In Ref. [34] , the electron neutrino-antineutrino ratio at source (Sec. 3.3) 
and detector (Sec. 4.3) has been computed explicitely. At the source, the 
following observations have been made: 

• Additional production modes (such as direct and multi-pion production) 
in photohadronic interactions produce 7r~ in addition to n + . The pure 7r + 
source therefore does not exist, and an up to 20% contamination from v e 
at the source has to be accepted even in this case. As a consequence, the 
Glashow resonance must be seen to some degree, even for the p^ source. 

• Since the photohadronic interaction in Eq. also exists for neutrons 
producing 7r~ in this case, any optical thickness to neutron escape will lead 
to a 7r + -7r _ equilibration. This means that the optically thick pj source 
cannot be distinguished from the pp source by the Glashow resonance. 

• Neutron decays, which are inherently present in any photohadronic source, 
lead to v e flux, faking a 7r~ contribution and therefore a pp source in 
particular energy intervals (determined by the maximal proton energy). 

As a consequence, only the "p7 optically thin" (to neutron escape) source 
might be uniquely identified if less than 20% of v e contamination are found. 
On the other hand, one cannot uniquely identify a pp source. Note that modern 
approaches take the composition between pp and pj interactions as a variable, 



see Refs. 89-91 . However, these approaches can typically not describe the 
contamination from neutron decays, because this depends very much on the 
model. 

At the detector, the use of the Glashow resonance becomes even more 
complicated for the following reasons: 

• Flavor mixing re-introduces a v e component from produced in /i + de- 
cays even for pure n + production, c/., Eq. 

• As a consequence, the flavor composition at the source is important, and 
has to be determined at the same time. For example, a muon damped pp 
source may be easily mixed up with a pion beam pj source. 

• The Glashow resonance occurs at a specific energy (6.3 PeV), which means 
that for this process the energy dependence is not important and only one 
particular energy matters. However, the transformation of the energy from 
source to detector depends on redshift and a possible Lorentz boost, c/., 



Eq. (15), which have to be known to draw any conclusions. 

In summary, the discovery of a neutrino signal with significantly suppressed 
Glashow resonance events will be interesting, and may allow for possible con- 
clusions about the source if flavor composition, z, and r are known. On the 
other hand, the detection of a pronounced Glashow resonance events is prob- 
ably non-conclusive for the physics of the source. 

4.4 Testing new physics in the neutrino propagation 



In the Standard Model, the transition probability P a p in Eq. (151 is described 



by the usual flavor mixing in Eq. (161, which is independent of energy. New 
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physics effects may lead to deviations from this picture, where the effects 
discussed in the literature include sterile neutrinos, neutrino decay, quantum 



decoherence, Lorentz invariance violation, among others, see Ref. 53 for a 



review. In some of these cases, P a p may even be a function of energy P a p(E). 
We choose neutrino decay as an example in this section, see Refs. 38 79 89 



92 - 991 , where in Refs. 97 - 99 energy dependent effects have been considered. 

Following Ref. 99 , for decay into invisible particletr^l the transition prob- 
ability can be described by a modified version of Eq. (|16|) : 



P, 



a/3 



3 

£ 

i=l 



\U m \ 2 D t (E) with Di(E) = exp 



A L 

-(Xi — 

l E 



(19) 



as the damping coefficient 



102 



Here = m^/r^ with rf the rest frame 



lifetime for mass eigenstate Vi. Typically the neutrino lifetime is emoted as 
rf /m.i since Wj is unknown, see Ref. 53 for a review. Note that this transition 
probability is e nerg y dependent, compared to Eq. ( 16 ). In this case, the flavor 
ratio R in Eq. (17) can be rewritten as 



R = 



P efl (E)X(E)+P^(E) 



[P ee (E) + P eT (E)} X(E) + [P^E) + P^(E)} 



(20) 



if X(E) = Q e (E)/Q ll (E) is the ratio between electron neutrinos and muon 
neutrinos ejected at the source (assuming that hardly any tau neutrinos are 



produced). Eq. (20) carries now two energy dependencies: the energy depen- 
dence of the flavor composition at the source X(E), and the energy depen- 
dence of the new physics effect P a /3(E), which have to be intrinsically disen- 
tangled. On the other hand, the energy dependence of the new physics effect 
may provide a unique signature |97||98| , and the energy dependence of the fla- 
vor composition at the source may help to disentangle different scenarios [99] . 

In general, there are 2 3 = 8 decay scenarios for invisible incomplete decays, 
since every (active) mass eigenstate may be either stable or unstable. In Fig. [9] 
left panel, R is shown as a function of the initial flavor composition X for these 
different (complete) decay scenarios, where filled disks correspond to stable 
mass eigenstates and unfilled disks to unstable mass eigenstates (the scenario 
with only unstable states is of course not shown, since no neutrinos can be 



See e.g. Refs. 



100 



for the systematical discussion/classification of possible 
scenarios. In general, neutrinos may decay into particles invisible for the detector, 
or other visible active flavors. Furthermore, the decay may be complete, i.e., all 
particles have decayed, or incomplete, where the energy dependence of the decay 
can be seen. For the case of complete decays, there are only eight different effective 
decay scenarios (each mass eigenstate can be stable or unstable), no matter if the 
decays are visible or invisible or if intermediate (unstable) invisible states are 
involved [89]. Likewise, there are only eight scenarios for invisible incomplete 
decays, whereas the treatment of incomplete visible decays is more complicated, 
see, e.g. Refs. [l0u][luT] . 
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Fig. 9. Left panel: R as a function of the initial flavor composition X for differ- 
ent decay scenarios (complete decays), where filled disks correspond to stable mass 
eigenstates and unfilled disks to unstable mass eigenstates. Right panel: R as a func- 
tion of energy for a pion beam to muon damped source (c/., Fig. [5j TP 13) for a 
specific decay parameter (same for all unstable states), as given in the panel. The 
dashed curves stand for complete decays. Figures adopted from Ref. [99|. 



detected in this case). In this panel, the different types of sources are also 
marked. One can clearly see that especially for a pion beam source, three of 
the scenarios cannot be distinguished, whereas these can be disentangled in 
principle for any other type of source. In addition, for three scenarios (the 
scenarios where only one mass eigenstate is stable), R is independent of the 
initial flavor composition. Note that the scenario with the largest R (mi and 
m2 unstable) faces the strongest constraints because of the observation of 
neutrinos from supernova 1987A. 

In Fig. (9} right panel, R is shown as a function of energy for a particular 
choice of the decay parameter times distance and a specific source with a flavor 
transition from pion beam to muon damped source at about 10 6 GeV. First 
of all note that because of the exponential damping in Eq. (19), the decays 
are practically complete for energies E <C 10 8 GeV, and the neutrinos are 
stable for energies E ^> 10 s GeV for the chosen decay parameters. This is why 
the scenarios deviate from the complete decay curves (dashed) at very high 
energies and converge into the standard scenario there (all neutrinos stable). 
The pion beam (low energies) cannot distinguish three of the scenarios at low 
energies, as expected. However, above 10 6 GeV, where the flavor transition 
occurs into a muon damped source, the corresponding curves split up (and 
the neutrino decays are still complete), before they converge into the standard 
case. This example illustrates that the energy dependent flavor transition of 
a specific source may be a useful for new physics tests, provided that enough 
statistics can be collected. 
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5 Application to generic (AGN-like) sources 

As it is illustrated in Fig. [6] (see also Eq. (|5j)), the proton N' p and photon densi- 
ties N' within the source control the secondary meson, and hence the neutrino 
production. In this section, we follow the ansatz in Ref. 34 : we assume that 



the target photons are produced in a self-consistent way, by the synchrotron 



emission of co-accelerated electrons (positrons); see, e.g., Ref. 103 for a cor- 
responding specific (BL Lac) AGN blazar model. The main purpose of this 
model is the prediction of spectral shape and flavor composition of the neu- 
trino source, while it cannot predict the normalization of the neutrino flux 
in the form presented here. In addition, no neutrinos have been observed yet, 
which means that an ansatz tailor-made for neutrinos may be useful to test the 
interplay between detector response and source model - after all, there may 
be sources for which the optical counterpart is absorbed, so-called "hidden 



sources", see, e.g., Refs. 104 107 . The ingredients to this model are compa- 
rable to the conventional GRB neutrino models, discussed in Sec. [6j where 
however the origin of the target photons is different. 



5.1 Additional model ingredients 

The primaries, in this case protons and electrons, are assumed to be injected 
with an (E')~ a injection spectrum, where the (universal) injection index a 
is one of the important model parameters. They are assumed to lose energy 
dominantly by synchrotron radiation, controlled by B' , and adiabatic cooling, 
controlled by R' , which also determine their maximal energies. This means 



that Eq. ( 13 ) is applied to the primaries, which implies that the acceleration 
and radiation zones are different in this model. As a consequence, the electron 
spectrum becomes loss-steepened by one power. 

The synchrotron photons (the target photons in Eq. ([5])) are computed 
in the Melrose-approximation |108| averaged over the pitch angle. The power 
radiated per photon energy e 1 by one particle with energy E' , mass m and 
charge q in a magnetic field B' is given by: 



W, *) = 1.8 ■ f" 5 ',, ■ (4) ' e -'/4 with < = »|» ( * 

sV ' ' 16e mch \e'J c 16m \mc 

We have to convolute this with the spectrum of radiating electrons as 



P>(e>) = / 
Jo 



dE' N' e {E')P' s (e',E'). (22) 



The number of produced photons per time can be computed with: 

P'(e') 

Q' i{ e') = t ^l. (23) 
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Fig. 10. Left panel: Possible acceleration sites in Hillas plot as a function of R 
and B (version adopted from M. Boratav). Right panel: Classification of sources for 
injection index a = 2 in this plot (see main text). Some points from left plot are 
shown for orientation, as well as two new points (12 and 13). Figure taken from 
Ref. [34] . 



Assuming that the photons escape (and hardly interact again), the steady 
photon spectrum, which is needed for the computation of photohadronic in- 
teractions, can be estimated by multiplying Q'^ with escape time t' CBC ~ R! j c 
of the photons. The synchrotron spectral index (s')~ a ~< obtained from this 
approach is a 7 ~ (a e — l)/2 + 1 = a/2 + 1, which means that the dependence 
on the primary injection index a is small. The rest of the computation follows 
Sec. H 

The main parameters in this model are a, R' , and B' . For the sake of 
simplicity, we assume in the following that the source is only moderately 
Lorentz boosted with respect to the observer's frame, i.e., R' ~ R and B' ~ B. 
In this case, a convenient description of the parameter space of interest is the 



Hillas plot 109 . In order to confine a particle in a magnetic field at the source, 
the Larmor radius has to be smaller than the extension of the acceleration 
region R. This can be translated into the Hillas condition for the maximal 
particle energy 

£max [GeV] ~ 0.03 -rj-Z-R [km] ■ B [G] . (24) 

Here Z is the charge (number of unit charges) of the accelerated particle, B 
is the magnetic field in Gauss, and r\ can be interpreted as an efficiency factor 
or linked to the characteristic velocity of the scattering centers. Potential 
cosmic ray sources are then often shown in a plot as a function of R and B, 



as it is illustrated in Fig. 10 left panel, by the numbered disks (see legend 



for possible source correspondences). Assuming that a source produces the 
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highest energetic cosmic rays with E ~ 10 20 eV, one can interpret Eq. p4[) 
as a necessary condition excluding the region below the r = 1 line in Fig. |10| 
(for protons with i] — 1). If one allows for relativistic boosts between source 
and observer, this condition is relaxed, as one can read off from the figure (in 
this case, R and B in the plot have to interpreted as R' and B'). However, 
this method does not take into account energy loss mechanisms, which lead 



to a qualitatively different picture, see, e.g., Refs. 110 111 , and which are 



implied in our model. In the following, we will study the complete parameter 
space covered by Fig. 10 (left) without any prejudice. Since the location of the 
sources in Fig. [lO] cannot be taken for granted, we will refer to the individual 
sources as "test points" (TP), and leave the actual interpretation to the reader. 

Concerning the limitations of the model, it certainly does not apply ex- 
actly to all types of sources. For example, in supernova remnants, pp (proton- 
proton) or pA (proton-nucleus) interactions may dominate the neutrino pro- 
duction, which would require additional parameters to describe the target 
protons or nucleons. In addition, at ultra-high energies, heavier nuclei may 
be accelerated. The spirit of this model is different: It is developed as the 
simplest (minimal) possibility including non-trivial magnetic field and flavor 
effects. Another ingredient is the target photon density, which is assumed 
to come from synchrotron emission of co-acceleratcd electrons here. In more 
realistic models, typically a combination of different radiation processes is 
at work. However, in many examples with strong magnetic fields, the spe- 
cific shape of the photon spectrum is less important for the neutrino spectral 
shape than the cooling and decay of the secondaries, which depend on particle 
physics only. To check this, we have tested the hypothesis that acceleration 
and radiation zones of the electrons are identical (i.e., the electron spectrum 
is not loss-steepened by synchrotron losses, which is actually a simpler version 
of this model) . Thus, while it is unlikely that the model applies exactly to a 
particular source, it may be used as a good starting hypothesis. 



5.2 Flavor composition at the source 

We discussed the flavor composition at the source already in Sec. |3.3| where 
we also showed several examples for the energy dependence in Fig. [5j Let us 
now approach this from a systematical point of view as a function of R and 
B (for a = 2). In this qualitative classification of the sources can be 

found in Fig. |10| right panel, where it is implied that a certain flavor ratio can 
be clearly identified over one order of magnitude in energy close enough to 



the spectral peak; see Ref. 34 . From this figure, we find that for B > 10 12 G, 
all charged species lose energy so rapidly that the neutron decays dominate 
the neutrino flux. For B > 10 10 G, several processes compete, leading to an 
undefined source. For B < 1 kG, the sources behave as classical pion beams, 
which typically applies to sources on galactic scales. In the intermediate range, 
1 kG < B < 10 10 G, the source classification somewhat depends on the 
spectral shape, since this affects possible muon pile- up effects, the energy range 
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close to the spectral peak, and the competition of several effects. Depending 
on R, muon beam to muon damped sources (as a function of energy), muon 
damped sources, and pion beam to muon damped sources are found. There 
is some dependence on a, which affects the spectral shapes. For instance, for 
a~4a similar pion spectrum as for GRBs is obtained, which leads to a pion 
beam to muon damped sources for typical GRB parameters in consistency 



with Refs. 37 38 . In summary, the pion beam source assumption is safe if 
B < 1 kG, whereas for stronger magnetic fields the magnetic field effects on 
the secondaries have to be taken into account. However, for given parameters, 
this flavor composition can be predicted. For instance, R may be estimated 
from the time variability of the source, B from energy equipartition arguments, 
and a from the observed photon spectrum. 

5.3 Interplay between spectral shape and detector response 



We discussed already earlier in Sec. |4.2| the interplay between spectral shape 
and detector response, see Ref. [46]. From Fig. [8] it is clear that spectral shapes 
with a peak at the position of the differential limit minimum can be better 
limited than others. This is especially clear if two event types or detectors are 
compared for the same spectrum. However, how should one quantify that for 
different spectra seen by the same detector? Consider, for example, the fluxes 
#2 and #11 in the left panel in Fig. [8] both leading to the same event rate by 
definition. Which of the two neutrino sources is the detector more sensitive to 
in terms of the physics of the source? In order to address this question, it is 
useful to assign a single number to each spectrum which measures how much 
energy in neutrinos can be tested for a specific spectrum and event type. We 
choose the energy flux density 

F l3 = J E v <f> p dE u (25) 

as this quantity, which we show in units of ergcm~ 2 s _1 for point sources in 
order to distinguish it from E^fi in units of GeVcm~ 2 s _1 (lerg ~ 624 GeV). 

This quantity measures the total energy flux in neutrinos, and it is useful 
as a performance indicator measuring the efficiency of neutrino production in 
the source. In order to sec that, consider the alternative derivation of Fr from 



Eq. (15) (neglecting a possible Lorentz boost of the source) 

F f> = E ^§r> where L -c=V I KQ'ujK (26) 



is the "neutrino luminosity" and V' is the volume of the interaction region. 
Since the neutrinos originate mostly from pion decays and take a certain frac- 
tion of the pion energy (about 1/4 per produced neutrino for each charged 
pion), the neutrino luminosity is directly proportional to the (internal) lumi- 
nosity of protons L[ nt (or the proton energy dissipated within a certain time 
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Fig. 11. Regions where the sensitivity in Fp exceeds 10 -9 erg cm -2 s _1 [sr _1 ] for 
P = n (left panel) and /3 = r (right panel) for several selected data samples (90% CL, 
a — 2). The dashed regions "UHECR" indicate where 10 20 eV cosmic ray protons 
are expected to be produced in the model. Figure taken from Ref. [46] . 



frame AT) and the fraction of the proton energy going into pion production, 
commonly denoted by f v (if the energy losses of the secondaries can be ne- 
glected). Since a possibly emitted photon flux can be often linked to L- lnt by 
energy partition arguments, one has Fp oc f„ x L mt oc f„ x L 7 , and Fp is a 
measure for f„ x L; nt of the source (if no photon counterpart is observed) , or 
even /„• itself (if a photon counterpart is observed). 

Regions for a specific sensitivity to Fp are shown as a function of R and B 



in Fig. 11 for (left panel) and v T (right panel), for several source declinations 
in IceCube and for Earth-skimming tau neutrinos in Auger (in this case, for 
a diffuse flux). There are several conclusions from this figure: first of all, 
IceCube reponds very well to the usual suspects, such as AGNs (left panel). 
Even for v T1 which produce a muon track in only 17% of all cases, most of 
this parameter space can be covered (right panel). Thus, it is clear that most 
sources will be also detected if the partition between and v T was heavily 
disturbed, such as by a new physics effect in the neutrino propagation. The 
sensitive region, however, somewhat depends on the source declination. For 
very large values of B, the neutrino energies are lower and instruments such 
as the DeepCore array may respond better to the flux (both panels). On 
the other hand, the region where the UHECR are expected in this model 
(lower right corner) is better covered by Auger (right panel) . This is not a big 
surprise: the neutrino spectrum follows the proton spectrum in the absence of 
strong magnetic field effects on the secondaries, which means that the spectra 
(c/., Fig. [8]) peak at high energies. A very interesting region may be the gap 
between the IceCube and Auger regions: perhaps a future instrument such as 
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KM3NcT should optimize their geometry to be complementary in terms of 
energy range coverage. 



6 GRB neutrinos and the multi-messenger connection 



Recall again that, as it is illustrated in Fig. [6] (c/., Eq. (5j)), the proton N' p 
and photon densities N' control neutrino production. In the previous section, 
we assumed that is produced by synchrotron emission of co-accelerated 
electrons. Here we emphasize the multi-messenger connection, assuming that 
N' can be derived from the gamma-ray observation 30||112 ; see also Refs. 



15 

the 



20 2l]|113 for analytical approaches. The motivation is alio different from 
previous section: whereas we were interested in a systematic parameter space 
study of spectral shape and flavor ratio before, the main emphasis here is the 
prediction of spectral shape, flavor compositon, and absolute neutrino flux 
normalization for a specific set of GRBs observed in gamma-rays. 

6.1 Additional model ingredients 

Here we describe the key ingredients of the conventional (numerical) fireball 
GRB model for neutrino emission, following Ref. [30], where we focus on the 
normalization. For a detailed comparison to the analytical calculations, see 
Refs. [29| 112] . Because the model does not describe the neutrino production 



in a time-resolved way, it makes sense to relate the neutrino production to 
the (bolometric) gamma-ray fluence Sboi (in units of ergcm~ 2 ) during the 
burst. The isotropic bolometric equivalent energy Ei SOj boi (in erg) in the source 
(engine) frame can be then obtained as 



4ndl 
(1 + *) 



S ho i . (27) 



It can easily be boosted into the SRF by E' iso bo] = Ei S0 ^ i/r. Assuming en- 
ergy equipartition between photons and electrons, the photons carry a fraction 
e e (fraction of energy in electrons) of the total energy -E'iso, tot and 

-E'iso, tot — e e 1 ' Eiso.bol • (28) 

In order to compute the photon and proton densities in the SRF, it turns out 
to be useful to define an "isotropic volume" V^ a 

v v ' s v ' 

.Rc Ad' 

where Rc is the collision radius and Ad' the shell thickness of the colliding 
shells. It can be estimated from the (observed) variability timescale t v , z, 
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and r, and can be regarded as the volume of the interaction region assuming 
isotropic emission by the source. Because of the intermittent nature of GRBs, 
the total fluence is assumed to be coming from N ~ Tgo/t v such interaction 
regions, where T90 is the duration of the burst (time during which 90% of the 
total energy is observed). Now one can determine the normalization of the 
photon spectrum in Eq. ^ from 



e'N'Je')de' 



EL 



,bol 



N V 

ISO 



(30) 



assuming that the spectral shape is determined by the observed spectrum. 
Similarly, one can compute the normalization of the proton spectrum in 
Eq. (§ by 

E'^E^AE'^^^^, (31) 



fe NV 



I 

ISO 



where f e is the ratio between energy in electrons and protons (/~ : baryonic 
loading). Note that in the end, we will obtain the neutrino flux <j> per time 
frame per interaction region from Eq. (151 with N — V^. Assuming that the 
magnetic field carries a fraction of E- lao tQt , one has in addition 
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Typical values used in the literature are f e < 
Ref. (2TI ). An explicit calculation for B' yields 
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(32) 



e e ~ (b — 0.1 (see, e.g., 



B' ~ 220 
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In summary, once the photon (from observation) and proton (typically 
E~ 2 ) spectral shapes are determined, the proton and photon densities in the 
source and B' can be calculated with the above formulas from the observables 
(gamma-ray fluence, r, z, t v ). Eq. (29) implies that for fixed t v , the larger 



r, the larger the interaction region, and the smaller the photon density in 



Eq. ( 30 ) , which directly enters the fraction of fireball proton energy lost into 
pion production f n cx -T~ 4 20 , or, consequently, E^(f> cx F~ 2 . Therefore, the 



main contribution to the neutrino flux is often believed to come from bursts 
with small Lorentz factors, see discussions in Refs. 20 114 115 . We follow 
this conventional fireball approach in the following. Note that this numerical 
approach contains all the ingredients of Refs. 20 21 explicitely - such as the 



cooling of the secondaries. The neutrino emission can be then easily computed 
as shown in Fig. [6] In addition, note that there are alternatives to the model. 
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For instance, if the bursts are alike in the comoving frame, as suggested in 



Ref. 116 , one has E 2 cj) oc F 2 [30] . One can read off from Eq. (29 1 that the 
correlation F oc t v is expected in that case, since V{ so will be similar for the 



bursts. See Ref. 30 for a discussion of the neutrino flux for different model 
hypotheses, and how the neutrino flux can in principle be used to discriminate 
among these. 

6.2 Systematics in the interpretation of aggregated fluxes 

Since the number of neutrinos expected from a single GRB is small, dedi- 
cated aggregation methods are needed. For instance, one may search for the 
diffuse flux from GRBs, which, however, has to fight the background from 
atmospheric neutrinos. Another possibility is to use the gamma-ray obser- 
vation to infer on time window and direction of the neutrino signal, which 
effectively leads to significantly reduced backgrounds. In addition, a proce- 
dure such as the one in Sec. |6.1| may be used to predict the absolute neutrino 
flux, its shape (see, e.g., Fig. |4j) , or its flavor composition (see, e.g., Fig. [7]) on 
a burst-by-burst basis. Summing over many observed bursts, such an analysis 



is also called stacking analysis, see, e.g., Ref. 15 for a recent example. The 
diffuse limit can be extrapolated from such a stacked flux, which is also called 
the quasi-diffuse limit; see Ref. [7] for details. Here we discuss some of the 
implications when a stacked neutrino flux is translated into a quasi-diffuse 
flux. 

In order to illustrate this problem, following Ref. [30], let us consider the 
redshift distribution of the GRBs as an example. We show in Fig. [12] left 
panel, a population of 10 000 bursts, representative for the number of GRBs 
in the visible universe over about 10 years (c/., histogram). These bursts are 



assumed to follow the star formation rate from Hopkins and Beacom 117 



with the correction £{z) from Kistler et al. |1 18 . For the sake of simplicity, 



assume that all bursts have the same isotropic luminosity. From the discussion 



after Eq. (15), we have E 2 (j> oc l/d^z) 2 , which means that closer GRBs, 
which are however rarer, will lead to a larger neutrino flux. Thus, it is the 
product d~[ 2 dN/dz which determines the main contribution to the neutrino 



flux, shown as solid curve in Fig. 12 While the peak contribution in terms 



of the GRB distribution is at z ~ 2 — 4, this contribution function peaks at 
z ~ 1. This observation has several implications: First of all, if the redshift z 
is not measured, the neutrino flux may be overestimated if z = 2 is assumed; 



c/., Eq. (27). Second, the number of bursts contributing in the region z ~ 1 is 
rather small, which means that large statistical fluctuations are expected in 
quasi-diffuse flux estimates based on the stacking of a few bursts only. 

We quantify this systematical errors in the right panel of Fig. [l2j In this 
panel, the probability that the quasi-diffuse flux extrapolated from a low 
statistics sample with n bursts is larger than a certain fraction of the dif- 
fuse limit (see legend for different values of n). This function corresponds to 
(one minus) the cumulative distribution function of the probability density, 
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Fig. 12. 



Left panel: Distribution of 10 000 bursts dN/dz as a function of redshift 

-2 



(histograms) and relative contribution of the individual GRBs d£ dN/dz (solid 
curves) . The dashed curve shows the exact distribution function. Here it is assumed 
that the GRBs follow the star formation rate from Hopkins and Beacom [117| with 
the correction £(z) from Kistler et al. 



118 . Right panel: Probability that the quasi- 



diffuse flux extrapolated from a low statistics sample with n bursts is larger than a 
certain fraction of the diffuse limit (see legend for different values of n). This function 
corresponds to (one minus) the cumulative distribution function of the probability 
density. Figure taken from Ref. [301. 



and the step function corresponds to the diffuse limit. One can read off from 
this figure that for n ~ 100, corresponding to the analysis in Ref. [15], the 
quasi-diffuse extrapolation will be within 50% of the diffuse limit in the prob- 
ability range corresponding to 90% of all cases (between 0.05 and 0.95). This 
means that a 50% error on the quasi-diffuse flux can be estimated from the 
redshift distribution only, while additional parameter variations increase this 
error 
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6.3 Neutrino flux predictions from gamma-ray observations 

Since IceCube has not observed any GRB neutrino flux yet, there has been 
increasing tension between the model predictions [9 20j 2l] and the obser- 
vation [l5j[l9](^] On the other hand, a direct comparison with the photo- 
meson production in Ref. [:9 was performed in Ref. p9] (see also Fig. [I]), 
demonstrating that the neutrino flux is actually underestimated in the ana- 
lytical approaches. Therefore, the differences between the numerical approach 
in Sec. 6.1 and the analytical models in Refs. 20 21] (based on Ref. (9|), have 



11 Recently, the superluminal propagation of neutrinos has also been proposed as a 
reason why no neutrinos have been seen, see, e.g.. 



Ref. 



124 . 
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Fig. 13. Expected (time-integrated) neutrino flux for three different GRBs (param- 
eters: see Refs. [Tl9l[720] for GRB 080916C, Refs. [TT9|[T2l] for GRB 090902B, and 
Refs. [Tl9p22] for GRB 091024). First row: Revision of IceCube analytical method 
CFB (correction of shape, normalization from pion production efficiency c/„, and 
normalization from neutrino versus proton spectral shape eg), leading to RFB. Sec- 
ond row: Comparison analytical (CFB, RFB) methods with simplified numerical 
method "WB Zi-approximation" Eq. (Jsj) and full P7 interactions. Third row: Impact 
of adiabatic cooling on protons and secondaries. Fourth and fifth row: correspond- 
ing inverse timescales (rates) for protons and secondaries, respectively. Courtesy of 
Svenja Hummer 11231. 
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been identified in Ref. 112 by a re-computation of the analytical models and 
the analytical computation of a simplified version of the numerical code. As 
far as the astrophysics ingredients are concerned, these approaches can be 
shown to be equivalent, based on the same logic; see Sec. |6.1| The main differ- 
ences are: magnetic field and flavor-dependent effects are explicitely included 
in numerical approach, additional pion, neutron, and kaon production modes 
are computed, and the full energy dependencies are taken into account. 

We illustrate in Fig. [13] the comparison between analytical and numerical 
approaches, where three different recent Fermi-measured GRBs have been 
chosen as examples: GRB 080916C, GRB 090902B, and GRB 091024. GRB 
080916C has been selected, because it is one of the brightest bursts ever seen, 
although at a large redshift, and one of the best studied Fermi-LAT bursts. 
The gamma-ray spectrum of GRB 090902B can be fit by a Band function and 
a cutoff power law (CPL), which means that it can be used to illustrate the 
difference. GRB 091024 can be regarded as a typical example representative 
for many Fermi-GBM bursts 119 , except from the long duration. Note that 
the first two bursts have an exceptionally large r > 1000, whereas r ~ 200 
for the third burst. All three bursts have in common that that the required 
parameters for the neutrino flux computation can be taken from the literature, 
in particular, the properties of the gamma-ray spectrum (including fluence), 
r, t v , z, and T90, see figure caption for the references. 

In Fig. 13 we show in the first row the computation of the predicted neu- 
trino flux with the IceCube analytical method 
fireball calculation) here. As described in Ref. 



21 , called CFB (conventional 



112 , the corrections of shape, 



normalization from pion production efficiency c f w , and normalization from 
neutrino versus proton spectral shape cs (see also Ref. [45]) lead to a revised 
analytical calculation RFB (revised fireball calculation). From the figure, one 
can easily read off that these revisions strongly depend on the burst param- 
eters, especially the photon spectral shape. Comparing the predicted CFB 
fluxes with the bursts used for the IC40 analysis (Fig. 1 in Ref. 15 ), one can 
easily see that the expected fluxes of the first two bursts are about a factor 
50 below that of the most luminous bursts in that analysis, and the third 
example about a factor of 5 below. This is expected from the scaling of the 
pion production efficiency cx .T~ 4 in that approach. 

In the second row of Fig. [13] we show the comparison of the analyt- 
ical (CFB, RFB) methods with a simplified numerical method "WB A- 
approximation" , c/., Eq. and the full ^7 interactions. In most cases, the 
simplified numerical approach matches the method RFB rather well, which 
proofs the validity of the derived corrections. However, for GRB 090902B 
(middle panel), the spectrum below the first break is different from the an- 
alytical estimate because the scalings of the weak decays limit the steepness 
of the spectrum there. The final numerical calculation including all produc- 
tion modes is then significantly enhanced again, especially due to multi-pion 
production [29]. Using a cutoff power law for GRB 090902B for the gamma- 
ray spectral fit (middle panel), the normalization of the prediction slightly 
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reduces in that example because the photon density above the photon break 
is suppressed. Comparing the original CFB method with the final numerical 
computation "full P7" , it is interesting that this can significantly deviate in 
both normalization and shape, but this deviation depends on the burst pa- 
rameters. Very interestingly, a similar neutrino flux normalization for all three 
bursts is obtained, which means that its probably not warranted to say that 
the neutrino flux from high-_T bursts such as GRB 080916C is expected to 
be small. However, note that such extreme bursts only make up for a small 
fraction of the observed bursts, and the conclusions from neutrinos will be 
determined by the statistical properties of the burst sample in the stacking 



analysis. Finally, note that the numerical calculations in Fig. 13 do not de- 
pend on any approximations, whereas different analytical methods lead to 
different predictions, similar to CFB. Therefore, the numerical computations 
should be regarded as the benchmark which defines the corrections, not vice 
versa. Within the simplest fireball neutrino model, there are only small model 
dependencies within the numerical approach. For instance, the integral limits 
in Eq. (31 1 (the minimal and maximal proton energies) have to be specified}^] 



and a bolometric correction may have to be applied to Eq. ( 30 ) - which typ- 
ically has small effects. 

One of the effects included in the final result of the second row of Fig. [13] 
is the adiabatic cooling of primaries (protons) and secondaries (muons, pi- 
ons, kaons). We illustrate the effect of this cooling component in the lower 
three rows of the figure, where we show the impact on the final result (third 
row) and the respective inverse timescales (rates) for protons (fourth row) 
and secondaries (fifth row). For the protons (fourth row), synchrotron losses 
are assumed to determine the maximal energies in the absence of adiabatic 
cooling, whereas the larger of the synchrotron or adiabatic cooling loss rates 
determines the maximal proton energy otherwise. From the fourth row, one 
can also easily read off that the energy losses due to p7 interactions are typ- 
ically sub-dominant. The comparison with the third row illustrates that, de- 
pending on the burst parameters, the kaon hump (the rightmost one) can be 
suppressed by the adiabatic cooling of the protons, whereas the normalization 
of the spectra is hardly affected. From the fifth row, one can also read off that 
adiabatic cooling may have a small effect on the muons, for which it some- 
times dominates in a small energy range. This leads to a small suppression of 
the first hump, coming from the muon decays (see left and middle panels of 
third row). 



We show in Fig. 14 the predicted quasi-diffuse neutrino flux from the above 
numerical method to the IC40 bursts for the same bursts and parameters (solid 
black curve) used in that analysis, which is about one order of magnitude lower 
than the IC40 limit and a factor of two below the current limit. In this figure 
a number of systematical errors are shown as well, such as the statistical 



12 Note, however, that for an (E' p ) 2 injection spectrum, the energy partition only 
logarithmically depends on the minimal and maximal proton energies. 
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Fig. 14. Numerical prediction (NFC) of the quasi-diffuse flux expected from the 
bursts used in the IC40-analysis, including the estimates for several model- or 
method-specific systematical uncertainties (see main text). In addition, the IC40 
limit is shown, and two expectations are shown for comparison (IC59+40 from 
and IC86 extrapolated for A 1 ^ 86 ~ 3 x A^ 40 from IC40; see, e.g., Ref. 
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Ref. 



Figure taken from Ref. 112 



error discussed in the previous subsection, and the estimated astrophysical 
uncertainty (by varying the unknown parameters, such as proton injection 
index a p = 1.8. ..2.2, variability timescale t v by one order of magnitude 
around the IceCube standard values, t v = 0.01 s for long bursts, r from 200 
to 1000, and the ratio e e /es from 0.1 to 10). In addition, note that z has only 
been measured for a few bursts used in the IC40 analysis, whereas z ~ 2 has 
been assumed for the long bursts with unknown z. As we illustrated in the 
previous subsection, this is potentially problematic, which means that a solid 
lower limit for the prediction can be only obtained for bursts with measured z 
(dashed black curve in Fig. 14). Note that our prediction varies not as strong 
as one may expect from /, a ij -1 r~ 4 20 . First of all, it is clear from Eq. ( 15 ) 



that E„(f> v (E) cx r because E v oc I '. Second, the synchrotron losses of the 



secondaries damp this variation 115 : for larger r, the energy densities in the 



source will decrease because of the energy equipartition, and consequently B' 
in Eq. (32 1. This reduces the energy losses of the secondaries, which means 
that more energy goes into the neutrinos. 

From the extrapolated IceCube limit for IC86 in Fig. [l4j it is obvious 
that IceCube will finally test the parameter space of the GRB fireball models, 
whereas the current limits already start to enter the meaningful parameter 
space. On the other hand, it is not clear what values of r and t v most bursts 
actually have, and what would contribute most to the neutrino flux. For ex- 
ample, a theoretical study of the dominant r to the diffuse neutrino flux has 
been performed in Ref. [30 for different model hypotheses. The conventional 
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fireball approach, which is presented above, leads to r ~ 200, which means 
that bursts such as GRB 091024 would dominate the neutrino flux. In this 
case, one would expect the prediction rather in the upper half of the shaded 
region of Fig. |14[ and a near-future detection in IceCube may be rather likely. 
On the other hand, recent studies for t„ seem to prefer larger average t v ~ 0.1 s 
than the t v ~ 0.01s used by IceCube for long bursts [125 , which points to- 
wards the lower half of the shaded region. Independent of these parameters, 
the neutrino flux prediction is proportional to l// e , which means that the 
final IceCube results will limit the baryonic loading as key parameter. 

Finally, note that similar considerations as for the GRBs may apply to 
AGN models, such as Refs. 10 12 , 126 , and it is yet to be seen what the 



impact on the cosmic ray connection is, see Ref. 22 



7 Summary and conclusions 

We have discussed the impact of particle physics effects on the neutrino fluxes 
from cosmic accelerators, based on a generic numerical approach of the pho- 
tohadronic charged meson production in sources optically thin to neutron es- 
cape; c/., Fig. [6] The starting point have been the proton and photon densities 
(spectra) within the source, which lead to the secondary meson production. We 
have included additional i-channel, higher resonance, and multi-pion produc- 
tion modes, neutrinos from kaon and neutron decays, the synchrotron cooling 
of the secondaries (pions, muons, and kaons), the helicity dependence of the 
muon decays, and the full spectral dependencies. While this approach may not 
be able to describe every source accurately, it can be regarded as the minimal 
approach to neutrino production by pj interactions using at least the well- 
known particle physics ingredients, to be extended by further components if 
applied to particular source classes. We have applied this framework to GRBs 
in Sec. |6.1[ where the target photons are inferred from the gamma-ray obser- 
vation. We have demonstrated that it contains all the ingredients in frequently 
used analytical models, such as Refs. [9^ [T5pO|[2T] , which means that it can be 
regarded as a numerical refined version of the conventional fireball neutrino 
model. In addition, we have discussed a generic AGN-like application where 
the target photons are produced by the synchrotron emission of co-accelerated 
electrons in Sec. [5] From the presentation in this review, it must be clear that 
this application is at a similar level as the GRB approach in terms of the 
ingredients. Applications to specific neutrino flux normalization predictions 
are, however, yet to be performed. 

The main consequences of the particle physics effects discussed in this 
review can be summarized as follows: 

• Magnetic field effects and flavor mixing change the flavor composition. The 
assumption of an v e : : v r flavor composition at the source of 1 : 2 : 
(pion beam) can only be justified for B' < 1 kG, while it will normally 
change as a function of energy for larger values of B' . 
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The neutrino spectral shape for p-y interactions depends on proton and 
photon injection spectra and the magnetic field and flavor effects. Roughly 
speaking, the E~ 2 assumption only holds for the special case {e')~ l for 
the photon spectrum. This means that this assumption is too simple for 
many sources, since the detector response depends on the interplay be- 
tween spectral shape and differential limit. For example, the detector may 
miss a neutrino signal because the spectrum does not peak at the right 
energy, although the fraction of energy going into pion production in the 
source is the same as in another detectable case. 

Additional pion production processes increase the neutrino production sig- 
nificantly, and lead to an impact on the spectral shape. Approximations 
using the Z\-resonance approximation are useful for analytical estimates, 
but they do not take into account the cross section dependence at high 
center-of-mass energies. The impact on the neutrino flux is a factor of a 
few. 

Spectral effects, such as the energy dependence of the proton interaction 
length or the energy dependence of the photon spectrum, lead a signifi- 
cant reduction of the GRB neutrino flux prediction of about one order of 
magnitude in comparison to analytical estimates. As a consequence, Ice- 
Cube has just entered the relevant part of the parameter space to test the 
simplest GRB fireball model. Similar effects on AGN flux models are yet 
to be tested. 

Neutrino fluxes from kaon and neutron decays are generically expected in 
any source. However, the neutrinos from neutron decays show typically 
up below the peak, the neutrinos from kaon decays above the peak. The 
relative importance of these effects depends on the strength of B' (which 
may separate the spectra from the different parents) and the maximal 
proton energy (which controls the neutron decay contribution). 
Using the Glashow resonance for the discrimination of pp and p7 inter- 
actions in the source is challenging. Especially, it is very hard to infer 
any physics information from the discovery of Glashow resonant events, 
whereas the absence of Glashow resonant events (in presence of a neutrino 
signal) can be interpreted as a pj source optically thin to neutron escape. 
Already from the current IceCube limits one can extrapolate that conclu- 
sions for astrophysical sources will mostly like be based on source classes, 
not individual sources. The aggregation of fluxes from many sources, how- 
ever, introduces new systematical errors. For instance, the extrapolation 
from the neutrino fluxes from 117 GRBs to a quasi-diffuse flux implies 
an error of at least 50% (at the 90% confidence level), from the redshift 
distribution only. 

In the presence of a high statistics neutrino flux (close to the current 
bound), useful information on possible new physics effects in the neutrino 
propagation may be inferred from the ratio between muon tracks and 
cascades. In this case, the energy dependence of the flavor ratio at the 
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source and the energy dependence of the new physics effects may provide 
the necessary information to identify the effect. 

In conclusion, in the presence of data, neutrino astrophysics is now at 
the point where rough analytical estimates for neutrino fluxes are not suffi- 
cient anymore, since in that case no reliable conclusions can be drawn for the 
astrophysical models and thus the origin of the cosmic rays. Especially the 
particle physics effects discussed in this review can be relatively easily taken 
into account, and they are well known. We have demonstrated with numer- 
ous examples that these effects on flux normalization, spectral shape, flavor 
composition, and neutrino-antincutrino composition cannot be neglected. For 
example, a correction of the GRB neutrino flux prediction of one order of 
magnitude has been identified. Therefore, a reliable treatment of the parti- 
cle physics of the neutrino production should be the state-of-the-art of any 
neutrino data interpretation in the multi-messenger context. 
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